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Abstract 



The critical properties of the 5 = 1/2 Heisenberg two- leg ladders are inves- 
tigated in a magnetic field. Combining the exact diagonalization method and 
the finite-size-scaling analysis based on conformal field theory, we calculate 
the critical exponents of spin correlation functions numerically. For a strong 
interchain coupling, magnetization dependence of the critical exponents shows 
characteristic behavior depending on the sign of the interchain coupling. We 
also calculate the critical exponents for the 5 = 1/2 Heisenberg two-leg ladder 
with a diagonal interaction, which is thought as a model Hamiltonian of the 
organic spin ladder compound Cu2(l,4-diazacycloheptane)2Cl4. Numerical 
results are compared with experimental results of temperature dependence of 
the NMR relaxation rate 1/T\. 
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I. INTRODUCTION 



The 5 = 1/2 bond alternating Heisenberg chains with the next-nearest-neighbor inter- 
action attract much attention. Schematic of this model is shown in Fig. |l](a). The alternat- 
ing nearest- neighbor interactions are controlled by Ji and J 2 , and the parameter J 3 is the 
strength of the next-nearest-neighbor interaction. Magnetic properties of the spin-Peierls 
compound CuGe03 is probably described by this model.0 Alternatively, as shown in Fig. 
|l](b), another aspect of this model is that this model describes a spin ladder system with a 
diagonal interaction; J 3 and J\ are strength of the interaction along the ladder and along 
the rung, respectively, and the diagonal interaction is represented by J 2 . It is thought that 
such a model describes the new spin ladder compound Cu2(l,4-diazacycloheptane)2Cl4 (ab- 
breviated CuHpCl) From the theoretical point of view, the plateau on the magnetization 
curve is observed in this model.ill 

This model is also related to the Haldane-gap system. For J 2 = 0, J 3 > and —J\j J3 S> 1 
or J 3 = 0, J2 > and —J1/J2 3> 1, two spins connected with the strong ferromagnetic 
interaction (Ji) form a triplet. Because it effectively works as a S = 1 spin, we can regard 
such a model as describing the S = 1 Heisenberg antiferromagnetic chain in the strong 
coupling limit .Bi 

Usually, such a system has a gap A between the singlet ground state and a triplet excited 
one. An external magnetic field lowers the energy of a triplet excited state, but does not 
change the singlet ground state energy. At a critical field H C1 (— A), the gap disappears and 
a transition from a non-magnetic state to a magnetic one occurs. In a stronger magnetic 
field H > H Cl , the system is gapless until the magnetization is saturated, unless the plateau 
appears on the magnetization curve. 

One of the important features of the organic compound is its small coupling constant. 
The inorganic spin ladder compound SrCu203, which is known for its excellent two-leg ladder 
character, has a gap A > 400K.i Thus, the critical field H Cl is not accessible. On the other 
hand, the organic compound CuHpCl which is thought to be the S = 1/2 Heisenberg two- leg 
ladder compound, has a smaller eap A ~ 1O.5K.0 Because the critical field H Cl ~ 7.5T can 
be accessed, some experiment si'i0~@ have been done in the magnetic state which has finite 
magnetization. 

On the basis of bosonization, Chitra and Giamarchii have studied such a gapped system, 
at the critical field H Cl and in its vicinity. They have concluded that for the S — 1/2 
Heisenberg two-leg ladder system, the critical exponent does not depend on a magnetic field 
in the region close to H C1 . Nevertheless, this is inconsistent with experimental data observed 
111 CuHpCl.0 Since CuHpCl has a strong interchain coupling, it is important to study a spin 
ladder system with a strong interchain coupling theoretically. 

In this paper, we will investigate the critical properties of the spin ladder system with or 
without a diagonal interaction (see Fig. [TJ(b)) in the magnetic state. In Sec. [H], we present 
a model Hamiltonian and the numerical methodBi developed by Sakai and Takahashi. 
Using this method, we calculate the critical exponent of the spin correlation function in the 
magnetic state for the S = 1/2 Heisenberg two-leg ladder with the strong antiferromagnetic 
interchain interaction in Sec. II|, and with the strong ferromagnetic one in Sec. fV[ Finally, 
we calculate the critical exponents taking a di ag onal interaction into account, and compare 
the numerical results with experimental datal!3 observed in CuHpCl in Sec. [V| The last 
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section is devoted to our conclusion. 



II. NUMERICAL METHOD 

We consider the following Hamiltonian, 

H = H Q + Hi, (1) 

N/2 N 

} + 2J 3 J2S l -S i+2l (2) 

i=i i=i 

N 

H 1 = -g» B H'£S?, (3) 

i=l 

where Si denotes the 5 = 1/2 operator in the z-th site, N is the total number of sites, and 
H is the strength of a magnetic field applied along z-axis. The distance between unit cells 
which consist of neighboring two sites is set equal to unity. Afterward we set J3 = 1 and 
g\iB = 1- The periodic boundary condition is applied. Since the system has a rotational 
symmetry about z-axis and a translational symmetry, we can classify the Hamiltonian into 
the subspace according to the magnetization M = J^S- and the wave vector k. The lowest 
energy in each subspace is calculated using the Lanczos algorithm. For the iV-site system, we 
define E^(N, M) as the lowest energy of Ho in the subspace specified by the magnetization 
M and the wave vector k. For given N and M, E^{N,M) takes the minimum at k — ko. 
From now on, we simply describe Ek (N, M) as E(N, M). Note that k = for even M and 
k = 7r for odd M, for the parameters used in this work. 

To studv critical properties of the system in a massless region, we apply the numeri- 
cal method0El developed by Sakai and Takahashi, which is summarized below. First, we 
consider the ground state energy. According to the conformal field theory for the one- 
dimensional quantum system,0 the ground state energy of the massless system depends on 
the system size as 

±E{N, M) ~ e{m) - A± (#-«>), (4) 

where e(m) is the ground state energy per site in the thermodynamic limit, m = M/N is 
the magnetization per site, v s is the sound velocity, and A = |cy s with c being the central 
charge. We approximate v s numerically to be, 

v ° = £Z IF ~ ]h~ko~\ {Ekl(N ' M) ~ E{N ' M)) ' (5) 

where k\ is the wave vector closest to ko, i.e. \ki — k \ = 4^. Using those values of A and 
v s , we can obtain the central charge c. 

Next, we consider some excitations to obtain the asymptotic form of the spin correlation 
function in the magnetic state. We define 5k as the difference of the wave vector between 
the ground state and the excited one. The spin excitation which increases the magnetization 
(M — > M + 1) has 5k = k c . In this paper, k c always equals it. Moreover, in the magnetic 
state, a gapless excitation can exist at the soft mode 5k = 2kp. Because the system is 
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gapless in the magnetic state, the spin correlation function should decay algebraically in the 
ground state. We define a spin operator 0; as a combination of the two S = 1/2 operators 
S21 and Soi + i. The asymptotic forms of the correlation function of the new spin operator 
should beO 

(0S<^> ~ cos(k c r)r-v (r-xx), (6) 
{4>o4> z r ) - Am 2 ~ cos(2k F r)r~ vZ (r -> 00). (7) 

According to the conformal field theory^! the critical exponents 77 and rf are obtained 
aS 

= E(N, M + 1) + g(jV, M — 1) — 2g(jV, M) 
/l ' E kl (N, M) — E(N, M) ' 1 ' 

E 2kF (N,M)-E(N,M) 
T,[IV E kl (N,M)-E(N,M)- [J) 

In the following sections, we calculate the critical exponents using Eqs. @ and (P). 



III. STRONG ANTIFERROMAGNETIC INTERCHAIN INTERACTIONS 

We study the S = 1/2 Heisenberg two- leg ladder system with a strong antiferromagnetic 
interchain interaction. The parameters we use are J\ = 5.0 and J2 = 0. At first, to examine 
whether Eq. (|J) is satisfied, we calculate size dependence of the ground state energy. In 
Fig. 0, we plot the ground state energy E(N, M)/N as a function of 1/N 2 for m — 1/8, 1/4 
and 3/8. Because these lines seem linear, we conclude that Eq. (f|) is satisfied. For each m, 
the ^-intercept of this line is the ground state energy in the thermodynamic limit and the 
gradient of this line is —A. From the values of A and v s , we can estimate the central charge; 
c = -A/v s . The estimated value of the central charge is equal to unity within a few percent. 
Therefore, the system belongs to the universality class of the Tomonaga-Luttinger liquid. 

We also calculate the critical exponents r\ and rf by Eqs. (||) and @. In Fig. || rj(N) 
and r) z (N) are shown as a function of the magnetization. Because r](N) and rj z (N) obviously 
depend on the system size, we plot ??(iV) and rj z (N) versus iV -2 for m = 1/4, 3/8 in Fig. |j. 
It seems that rj(N) and rj z (N) are proportional to iV -2 except rf for m = 1/4. We evaluate 
the values of rj and rf in the thermodynamic limit as the y-intercept of the lines obtained 
by the least-squares method, except vf for m = 1/4. The deviation from the line is so small 
(less than 0.1%) that we neglect it. For vf for m = 1/4, we extrapolate using the only two 
values for the largest system and for the next-largest one. For other m, although rj(N) and 
rj z (N) can be calculated only for two or three values of N, we presume that rj(N) and i] z (N) 
are proportional to iV -2 . In Fig. |5], rj and rf in the thermodynamic limit are shown. 

If the system is described by the Tomonaga-Luttinger liquid, the universal relation i]-r] z = 
1 must be satisfied.Hl The values of rj ■ r] z in the thermodynamic limit are shown in Table 
|. From these values, we conclude that the relation rj ■ rf = 1 is satisfied well. This is 
consistent with c = 1 obtained above. We note that the value of rj ■ rf for m = 1/4 is 
slightly smaller than that for other m. This may be due to the peculiar size dependence of 
r) z (N) for m = 1/4. At m = 1/4, the plateau exists on the magnetization curve when the 
coupling constants are selected adequately! This means the system has a gap and Eq. (|4]) is 
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not satisfied. Nevertheless, for the parameter we use here, there must not be the plateau on 
the magnetization curve at m = 1/4. To declare the relation between the size dependence 
of r] z for m = 1/4 and the plateau on the magnetization curve, further studies are needed. 

IV. STRONG FERROMAGNETIC INTERCHAIN INTERACTIONS 

Next, we consider critical exponents for the spin ladder with a strong ferromagnetic 
interchain interaction, using J x = —5.0 and J 2 = 0. This system is expected to behave like 
the Haldane-gap system! because the strong interchain coupling makes two spins triplet 
along a rung and it works as one S = 1 spin effectively. Since we cannot obtain the value of 
E2k F accurately, we calculate only rj in the magnetic state. The magnetization dependence 
of r) is shown in Fig. || Because the size dependence of rj is small, we neglect the size 
correction. This result is in good agreement with the value of r? for S = 1 antiferromagnetic 

173 to to / to 

chainEJ not only qualitatively but quantitatively. 

As mentioned above, we cannot calculate rf directly. However, if the system is described 
by the Tomonaga-Luttinger liquid and the central charge c is equal to unity, r\ ■ rf = 1 
must be satisfied also in this system. Thus, we can estimate rf = r]^ 1 . To ascertain 
c = 1, we calculate size dependence of the central charge and extrapolate the result to the 
thermodynamic limit. For the iV-site system, we define A(N) and v s (N) as follows. Let 
N' be the greatest integer less than N so that m = M/N = M'/N'. E(N,M)/N and 
E(N' , M')/N' are calculated numerically. Fitting these two values with Eq. (|4]), we obtain 
A(N) as the gradient of this line. We also define v s (N) as v s for finite N, using Eq. (^j). 
Consequently, the central charge of the iV-site system is defined as c(N) = -A(N)/v s (N). 
For m = 1/4 and 3/8, the size dependence of c(N) is shown in Fig. [7[ Since c(N) seems to 
be proportional to N~ 2 , we fit these data with a straight line using the least-squares method 
and define c as the y-intercept of this line. The results are c = 0.995 for m = 1/4 and 
c = 0.998 for m = 3/8. These values suggest that c = 1. For the S = 1 antiferromagnetic 
chain with finite magnetization (0 < m < 1), the central charge is unity and the relation 
rj ■ r] z = 1 is satisfied.E-3 Since the value of rj obtained here is in good agreement with that 
for the S — 1 antiferromagnetic chain, the magnetization dependence of if in this system 
must have similar behavior to that for the 5 = 1 antiferromagnetic chain. 

Comparing these two results (Fig. [| and Fig. |6|), magnetization dependence of critical 
exponents shows characteristic behavior depending on the sign of the interchain interaction 
J\. Similar behavior was found by Sakai0 for the 5 = 1/2 bond-alternating chains, which 
is described by the Hamiltonian ([]]) restricting J 3 to zero and J\ to unity. He has concluded 
that this system is described by the Tomonaga-Luttinger liquid and its critical exponents 
satisfy rj z /2 < 1 < 2i] for J2 > and 2r] < 1 < r] z /2 for J2 < 0, in the magnetic state. 
Furthermore, 2r] = r] z /2 = 1 in the limits of m — * 0+ and m — > 1/2—, regardless of the sign 
of J 2 . 

V. SPIN LADDER MATERIAL: Cu 2 (l,4-diazacycloheptane) 2 Cl 4 

As mentioned previously, CuHpCl is described well by the isolated coupled S = 1/2 
chain so far. Recently, Chaboussant et al.0 have presented NMR study of this material 
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in magnetic fields. They have measured the proton NMR relaxation rate 1/Ti, varying 
magnetic fields. In a magnetic state, they have observed the divergence of 1/T\ at low 
temperature. If the system is gapless, 1/T\ diverges algebraically at low temperature, 

1 oc T~ a . (10) 
J-x 

They have concluded that a equals 1/2 at m — > 0+ and increases with magnetic fields till 
certain strength of magnetic fields. After that, an increase in the magnetic field weakens the 
divergence. Anyway, a increases with the magnetic field near H Cl (for small m). As they 
have pointed out, the exponent a is related to the exponent r] as rj = 1 — a. Hence, rj must 
decrease with an increase in the magnetic field for small m. 

It is believed that this material is a model system in the strong coupling (J1/J3 ~ 
Thus, the behavior of rj must be similar to that of the spin ladder system with a 
strong antiferromagnetic interchain interaction as shown in Sec. [TTT| . However, in Fig. the 
exponent rj increases with an increase in H from H Cl . This means that our numerical results 
differ from the experimental findings. 

It is thought that CuHpCl has not only ordinary ladder interactions J\ and J3 but a 
diagonal interaction J 2 .&^ To decide coupling constants, Hay ward et al.S have calculated 
magnetization curves and fitted them on the experimental data. They have suggested that 
the diagonal interaction is ferromagnetic and its value is J2/J1 = —0.1 for J3/J1 = 0.18. 
Hence, we calculate rj and r] z again, turning J2 on. The parameters we use are J\ = 5.5 and 
J2 = —0.55, corresponding to the above ratio. The values of rj(N) and rj z (N) obtained from 
Eqs. (|8]) and (Q) have size dependence. Thus, we extrapolate to the thermodynamic limit 
by the same way in Sec. |T|. Note that all of the rj(N) and rj z (N), including rj z for m = 1/4, 
are proportional to N~~ 2 when m is fixed. The numerical results are shown in Fig. |S] and the 
values of 77 ■ rj z are shown in Table |I[ Also in this system, the relation rj ■ rj z = 1 is satisfied 
well. 

In Fig. M, the maximum (minimum) of i] (r] z ) is smaller (larger) than the results in 



Sec. |HI[ Further, the value of magnetization at which rj takes the maximum shifts to a 
larger value region. Nevertheless, the ferromagnetic diagonal interaction does not change 
the magnetization dependence of rj and 77^, qualitatively. Our calculation still differs from 
experimental results. 

On actual materials, more complicated mechanisms may exist, e.g. interladder cou- 
pling, lattice distortion and so on. Recently, Nagaosa and MurakamiB3 have treated the 
S — 1/2 Heisenberg two-leg ladder system taking a lattice distortion into account. They 
have concluded that the lattice will be distorted to modulate the interchain coupling with 
the incommensurate wave vector which is proportional to the magnetization. To explain ex- 
perimental results, such effects may be important. Moreover, the coupling constants are still 
controversial.0 To determine the coupling constants, the value of magnetization at which rj 
takes the maximum may be useful, because the experimental determination of the absolute 
value of the critical exponent is formidable. 
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VI. CONCLUSION 



We have calculated the critical exponents of the spin correlation functions of the S — 1/2 
Heisenberg two-leg ladder system with a diagonal interaction in a magnetic field. First, 
we have considered pure spin ladders (a diagonal interaction is set to zero). Depending 
on the sign of the interchain coupling, the magnetization dependence of the critical expo- 
nents exhibits characteristic behavior. For a ladder with a strong ferromagnetic interchain 
interaction, the critical exponents depend on the magnetization like that for the 5 = 1 
antiferromagnetic chains not only qualitatively but quantitatively. Next, we have taken a 
ferromagnetic diagonal interaction into account in connection with CuHpCl. This diago- 
nal interaction makes the maximum of 77 small and the minimum of i] z large. Also, the 
magnetization at which 77 (r] z ) takes the maximum (minimum) becomes large. Comparing 
these results with experimental data, there is qualitative difference. To explain experimental 
results, more complicated effects may be important. 
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FIGURES 



FIG. 1. Schematics of the model considered in this paper. 

FIG. 2. Plots of E(N,M)/N vs. 1/N 2 with m fixed for 1/8, 1/4 and 3/8. The origin is 
shifted along the vertical axis. The values of points A and B are as follows: A = —2.82, —1.56, 
0.0; B = -2.77, -1.51, 0.05 for m = 1/8, 1/4 and 3/8, respectively. The dashed lines are to guide 
the reader's eye. 

FIG. 3. Magnetization dependence of the critical exponents r\ and rf for various TV. The 
parameters used here are J\ =5.0 and J 2 = 0. The dashed lines are to guide the reader's eye. 

FIG. 4. Size dependence of 77 (upper two figures) and rf (lower two figures) with m fixed. The 
horizontal axis means iV~ 2 . The dashed lines are obtained using the least-squares method and the 
dotted line is to guide reader's eyes. 

FIG. 5. Magnetization dependence of the critical exponents 77 and rf for the S =1/2 two-leg 
ladder with a strong antiferromagnetic interchain interaction ( J\ = 5.0, J2 = 0) in the thermody- 
namic limit. The dashed lines are to guide the reader's eye. 

FIG. 6. Magnetization dependence of the critical exponent rj for the 5 = 1/2 two-leg ladder 
with a strong ferromagnetic interchain interaction (Ji = — 5.0, J2 = 0). The dashed lines are to 
guide the reader's eye. 

FIG. 7. Size dependence of the central charge c. The dashed lines are obtained using the 
least-squares method. 

FIG. 8. Magnetization dependence of the critical exponents 77 and rf for the S = 1/2 two-leg 
ladder with a diagonal interaction ( J\ = 5.5, J2 = —0.55). The results for some N and extrapolated 
values are shown. The dashed and solid lines are to guide the reader's eye. 
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TABLES 



TABLE I. The values of r\ • rf in the thermodynamic limit for the 5 = 1/2 two- leg ladder with 
a strong antiferromagnetic interchain interaction (Ji = 5.0, J2 = 0). 



m 


1/12 


1/8 


1/6 


1/4 


1/3 


3/8 


2/5 


5/12 


7/16 



.993 .998 .997 .977 .999 .999 .999 .999 1.000 



TABLE II. The values of r] • rf in the thermodynamic limit for the S = 1/2 two- leg ladder 
with a diagonal interaction (J\ = 5.5, J2 = —0.55). 



m 


1/12 


1/8 


1/6 


1/4 


1/3 


3/8 


2/5 


5/12 


7/16 


rj ■ rf 


.996 


.999 


.999 


.995 


1.000 


1.000 


1.000 


1.000 


1.000 
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